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ABSTRACT 

A  mechanism  Is  proposed  as  a  means  of  passively  leveling  an  optical 
assembly  about  two  horizontal  axes  within  a  tolerance  of  one  arc-second. 

,  The  mechanism  Is  designed  so  that  It  will  stimulate  a  spherical  pendulum 

for  Input  motions  up  to  100  arc-seconds.  In  order  to  Insure  the  operation 
of  the  mechanism  and  establish  the  proper  geometry,  three  analytical 
models  were  developed.  Planar  and  spatial  models  which  assumed  rigid 
links  with  pin  Joints  and  a  planar  model  assuming  deformable  links  with 
pin  and  clamped  Joints  were  studied.  These  results  are  compared  with 
|  experimental  data  gathered  for  three  geometry  combinations. 

It  Is  shown  that  the  results  from  the  planar  and  spatial  analytics 
agree  to  within  10%,  with  the  spatial  model  predicting  less  error.  Hence, 
j  the  less  complicated  planar  analysis  may  be  performed  with  little  loss  In 

accuracy.  Neither  rigid  link  models  compare  well  with  the  experimental 
data,  with  as  much  as  an  80%  discrepancy  In  some  cases.  When 
j  deformable  links  are  Included  In  the  planar  model,  the  results  are  within 

16%  of  the  experimental  results.  This  Indicates  that  the  resistance  to 
bending  at  the  points  can  not  be  neglected  In  modeling  the  mechanism. 
|  The  mechanism  to  act  like  a  spherical  pendulum,  and  guidelines  are 

established  for  the  design. 

i 

I 
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Fig  1.1a  daplcta  the  elements  of  tuts  coactpt.  The 
pallet  earriaa  a  4  x  4  x  4  matrix  of  detectors  which  are 
targatad  by  a  laaar.  An  antira  row  of  datactora  la  targeted 
slnultaneoualy  by  a  borlsontal  tbln  plana  of  light.  Tbla 
plana  of  llgbt  la  created  by  aendlng  tba  laser  bean  through 
a  vertically  nountad  cylindrical  lana  of  abort  focal  length 
aa  In  Figure  1.1b.  An  entire  colunn  of  detectors  la  tar* 
gated  simultaneously  by  rotating  the  cylindrical  lens  Into  a 
borlsontal  position  as  In  Fig  1.1c. 

Figures  1.2  and  1.3  show  the  emitter  system  In  more 
detail  and  Its  passive  leveling  suspension  which  Is  Intended 
to  hold  the  Instrument  level  to  within  1  arc  sec. 

1 .2  Variables  In  a  Precision  Two  Axis  Level 

Figure  1.2  Is  an  Idealised  suspension  of  flexible  file* 
wants  which  converge  to  a  point.  Such  a  pendulum  will  pro¬ 
duce  a  level  platform  If  the  C.G.  Is  directly  below  the  ver¬ 
tex  of  the  system  and  If  all  the  suspension  elements  are 
equal  length  and  have  Infinite  flexibility.  In  reality, 
none  of  these  conditions  Is  met.  In  the  light  of  the  error 

goal  o£  8  *  1  arc  second,  the  minute  deviations  from  the 
0 

Ideal  geometry  can  have  significant  effects.  To  begin  with, 
there  Is  no  true  vertex  to  the  suspension  filaments  because 
the  wires  must  be  at  least  .006  In  diameter  to  support  even 
the  low  weight  of  the  lasers  and  lenses.  The  Inset  of  Fig¬ 
ure  1.2  demonstrates  that  Instead  of  a  vertex,  the  wires  are 


Actually  asperated  by  a  distance  6  ao  that  the  auapanalon  is 

of  the  kinematic  type  shown  by  Fig  1.3.  The  upper  plate 

represents  the  housing- of  the  instrument.  When  the  housing 

is  tipped  through  an  angle  9^  the  bottom  platform  will  not 

remain  level  but  will  experience  a  tilt  8  .  Although  8  is 

0*0 

smaller  than  8^,  it  could  nevertheless  exceed  the  tight 
objective  of  8q  <  1  arc  sac.  The  disruptive  role  of  unequal 
wire  lengths  and  tolerances  in  C.G.  location  will  contribute 
to  the  errors  in  the  instrument.  The  purpose  of  the  initial 
research  was  to  datarnlna  the  role  of  these  tolerances  to 
the  errors  in  instrument  suspensions. 


It  is  doubtful  that  the  instrument  will  be  mounted  on  a 

stationary  surface  since  it  will  be  used  in  a  production 

environment.  It  is  therefore  essential  to  minimize  any 

oscillations  of  the  optics  in  order  to  maintain  the  tight 

constraints  previously  Imposed  on  the  system.  So  as  not  to 

interfere  with  the  mechanism,  the  damper  must  be  free  of 

friction.  Figure  1.4  illustrates  an  air  damper  consisting 

of  two  cylinders,  the  inner  attached  to  the  bottom  plate  of 

the  mechanism  while  the  outer  cylinder  is  attached  to 

ground.  Although  simple  In  design,  oscillations  in  both  the 

X  and  Z  directions  can  be  effectively  damped.  Air  damping 

is  more  difficult  to  achieve  due  to  the  low  viscosity  of  the 

air,  while  on  the  other  hand,  it  does  not  require  seals  or 

* 

submersion  of  the  optics  as  would  be  required  with  oil  or 


.Spherical  Joints 
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Figure  1.3 


Schematic  representation  of  the  instrument 
indicating  a  3  dimensional  mechanism. 


Cylinder  Attached 
to  Link  5 


x  X 


Cylinder  Attached 
to  Ground 


Figure  1.4  Cylindrical  air  damper  for  the  two  axes  level 


The  last  section  of  this  report  focuses  on  the  snslysis 
end  expert me ntstlon  required  to  specify  sppropriste  geometr¬ 
ies!  rels tionships  to  fscilitste  the  design  of  such  a 


dsmper . 


CHAPTER  2 


PLANAR  ANALYSIS 


2.1  Obj  tc tlvii 

Tb«  goal  of  tills  analysis  and  those  to  follow  Is  to 

datarnlne  the  Influence  of  ttia  geometry  of  the  mechanism  on 

positioning  errors.  The  positioning  error  is  defined  to  be 

a  deviation  of  the  bottom  plate  from  the  horizontal.  This 

angle  will  be  represented  throughout  as  0  .  The  results  of 

o 

this  analysis  will  be  used  to  determine  the  feasibility  of 

this  type  of  mechanism  to  remain  level  and  the  geometric 

constraints  which  must  be  satisfied  in  order  to  minimize  0  . 

o 

2.2  Development  of  Equations 

The  mechanism  under  investigation  is  a  3  dimensional 
linkage.  However,  the  problem  will  be  addressed  initially  by 
reducing  the  generel  case  to  2  dimensions.  The  simpler  case 
offers  expedient  insight  into  the  significant  parameters  of 
the  system,  which  for  this  analysis  will  be  the  lengths  of 
the  input  link  (6)  and  output  link  (L),  height  of  the 
mechanism  (h),  the  centering  of  the  optics  (A)  and  the  die* 
tance  that  the  optics  lie  in  a  vertical  plane  from  the 
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output  link  (e).  The  2  dimensional  mechanism  is  shown  in 


Figure  2.1. 

The  sssunptlons  used  in  this  enalysis  are: 

*  The  cables  are  lnflnitly  flexible 

*  Frictionless  pin  joints 

*  Constant  link  lengths 

*  Links  2  &  4  are  two  force  nesbars. 


The  modal  Is  reduced  to  a  four-bar  linkage  with  tne 
Input  rotation  being  link  1,  as  shown  In  Figure  2.2.  A  mass 
Is  fixed  to  link  3  at  soma  known  location,  G,  relative  to 
pin  joint  A.  The  kinematics  of  the  linkage  can  be  written 


as . 

(*AC>2  -  (Ax-Cx)2  +  (Ay-Cy)2  (2.1) 

UBD>2  -  (Bx-Dx)2  +  (By-Dy)2  (2.2) 

UDC)2  -  (Dx-Cx)2  +  (Dy-Cy)2  (2.3) 

The  location  of  the  center  of  gravity  is  fixed  by  link  5. 

(Agc)2  -  (Gx-Cx)2  +  (Gy-Cy)2  (2.4) 

Ugd>2  "  <°*-Dx>2  +  (Cy-Dy)2  (2.5) 


Referring  to  Figure  2.2,  for  static  equilibrium  the 
forces  and  moments  on  link  3  are  conserved  which  yields: 


E  F  -  0 
3 


'acV/c  +  'bd'./b  -  *‘J  ■  0 


<2.0 


Where , 


E  »c  -  0 

TD/C  *  *»d7»/d  '  *C/C  x  “J  ’  0 


~  .  (Ax-Cx>y  +  (Ay- Cy )— 

A/C  J 


AC 


AC 


7  „  ( Bx-Dx)—  (Bjr-Py)- 

»'»  *BD  *BD 

7  ,  .  <c«ic«>i  +  (GjL-.cy>7 

G/C  1  1  J 

GC  GC 


Insertion  of  the  unit  vectors  Into  Equation  2 


-p^(Ax-Cx)  +  -p£(Bx-Dx)  -  0  si 
XAC  nD 


X  A  f*  X  BD  — 

-j - (Ay-Cy)  +  ** — (By-Dy)  -  mg  •  0  :j 

AC  *BD 


and  Equation  2.7  becomes 


j(Dx-Cx)(By-Dy)-(Dy-Cy)(Bx-DxA  -mg(Gx-Cx)  -  0 
XBD  1  ^ 


Equations  2.L  -  2.5  and  2.8  -  2.10  describe 
unknowns  for  the  position  solution  of  link  3 
points  A,  B  and  the  geometry,  links  1,  2,  3  and 
solution  was  carried  out  by  solving  the  eight 
equations  for  points  C,D  and  G  numerically  using 


12 

(2.7) 


.6  yields: 

(2.8) 

(2.9) 

(2.10) 

the  eight 
(given  the 
4).  The 
non- 1  inear 
a  Newton- 


Rapson  Method 


The  ouput  angle  9q  can  then  be  found  from  the  tangent  of 
points  D  and  C. 


2.3  Results  from  Planar  Analysis 


In  order  to  generalize  the  results,  the  geometry  will 

be  described  as  ratios  of  link  3  or  "L".  This  analysis 

showed  that  for  small  angles,  the  input  (6  )  output  (6  ) 

i  o 

relationship  was  linear  as  would  be  expected  for  small 

e 

O 

motions.  The  ratio  of  rr—  has  been  labeled  K  and  is  tfte 
measure  of  the  performance  of  the  mechanism. 


The  K  ratio  is  plotted  against  the  ~  ratio  in  Figure 
2.3  for  several  cases  of  £■.  As  can  be  seen,  K  goes  to  zero 

la 

only  when  the  input  link,  (6),  goes  to  zero  (a  true  pendu¬ 
lum)  or  the  height  h  goes  to  infinity  (the  effect  of  6  is 
lessened).  But  there  is  a  physical  limitation  of  the  size 
of  both  6  and  h.  The  construction  of  the  top  link  forces  6 
to  be  at  least  one  wire  diameter.  (See  Insert  on  Figure 
1.2)  Therefore,  the  only  way  to  decrease  ~  is  to  increase 
the  length  of  the  output  link.  For  packaging  purposes,  It 
is  desirable  to  minimize  the  height  of  the  mechanism.  From 

Figure  2.3,  it  is  therefore  advisable  to  generate  a  mechan- 

h  n 

ism  which  has  an  —  ratio  about  unity  since  —  much  greater 

L  la 

than  unity  offers  little  improvement  on  K. 


The  effect  of  Is  illustrated  in  Figure  2.4,  where  it 
A 

can  be  seen  that  —  may  vary  up  50%  without  introducing  posl- 
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Figure  2. 
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K  vs  6/L  es  e  function  of  h/L 
from  the  planar  analysis. 
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tioning  errors.  At  first  this  result  seemed  spurious,  but  e 

look  et  the  couf igure tion  of  the  mecnealsm  with  9  ■  0 

o 

explains  the  results.  The  inset  of  Figure  2.4  shows  that 
during  asseably  the  instant  center  of  link  1  with  respect  to 
link  3  is  brought  to  a  point  directly  above  the  ness  center 
of  the  nechanism,  while  the  wire  lengths  are  altered  such 
that  the  bottom  link  is  horizontal.  By  a  simple  tri¬ 
gonometric  analysis,  it  can  be  shown  that  this  assembly  pro¬ 
cedure  causes  the  distance  from  the  center  of  mass  to  the 
instant  center  to  be  independent  of  A.  (The  distance 
between  links  1  and  3  is  fixed  by  h.) 

Since  the  angular  velocity  ratio  tor  any  two  links  in  a 
mechanism  is  inversely  proportional  to  the  distances  from 
the  Instant  centers  in  the  frame  about  which  the  links  are 
rotating  to  the  instant  center  which  is  common  to  the  two 
links,  the  assembly  procedure  fixes  the  velocity  ratio.  The 
angular-  velocity  ratio  is  simply  K  since  the  relative  rota¬ 
tions  are  small  in  this  mechanism.  Thus,  K  remains  unal¬ 
tered  because  the  distances  from  the  Instant  center  to  links 
1  and  3  are  independent  of  A  and  therefore  fixing  the  velo¬ 
city  ratio,  K,  by  this  constant  distance  ratio. 

e 

The  effect  of  —  on  K  is  shown  in  Figure  2.5.  It  is 
Important  to  note  that  the  mass  center  should  lie  below  the 


bottom  link,  since  K  increases  as  the  output  link  length  e 
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decreases.  With  *  negative  value  of  e, 

lcally,  whereas  a  large  positive  value  f 

e 

reduces  K.  Since  the  value  of  —  Is  a 
recommended  to  make  this  value  as  large 


K  Increases  monoton- 
or  e  significantly 
design  choice,  It  Is 
as  possible. 


to  b«  of  any  value.  The  vector  loop  approach  products  cou¬ 
pled,  transcendental  equatlona.  While  this  creates  no  diffi¬ 
culties  in  planar  analysis,  standard  nusserical  techniques 
experience  difficulties  in  converging  on  a  solution  to  the 
larg.  nuaber  of  equations  required  for  spatial  analysis. 
The  situation  does  not  change  when  saall  rotations  are 
assumed  about  the  equilibrium  position  because  the  equations 
remain  linearly  coupled  by  cross  products.  Since  the  varl- 


11,  again  convergence  problems  arise  in  numeri¬ 


cal  equation  solvers. 


Minimisation  of  potential  anergy  was  explored  as  an 
alternative  means  of  developing  the  equations  of  motion. 
The  creation  of  multiple  partial  differentials  which  could 
not  be  determined  algebraically  forced  the  use  of  Lagrange 
Multipliers  to  evaluate  the  partlals.  When  the  Lagrange 
Multipliers  are  coupled  with  the  constraint  equations,  the 
number  of  equations  is  Increased  to  30. 


The  final  approach  contained  herein  is  less  elegant  in 
nature  out  requires  significantly  fewer  constrelnt  equations 
(13).  These  equations  are  still  non-linear,  but  the  equa¬ 
tion  solver  converges  rapidly  because  the  variables,  the 
coordinates  of  the  spherical  joints,  are  relatively  large. 


Mechanisms  Program),  A  Computer-Aided  Design  Analysis 
System  for  Mechanisms  and  Linkages",  Journal  Engineering 
for  Industry,  Trans.  ASME,  Ser.  B,  Vol.  99,  no.  3, 
August  1977,  pp.  773-778. 


Tills  approach  to  the  spatial  analysis  is  basically  the 
sane  as  the  planar  analysis,  except  it  allows  motion  in 
three  directions;  that  is,  the  mass  of  the  optics  is  at  some 
point  G  on  the  ternary  output  link  (as  shown  in  Figure  3.1) 
and  the  link  lengths  are  constant  with  frictionless  spheri¬ 
cal  joints.  The  input  rotation  is  about  the  X  and  Z  axis 
located  at  the  center  of  the  input  link.  Rotation  of  the 
input  link  about  the  Y  axis  does  not  alter  the  magnitude  of 
the  output,  in  that  it  merely  rotates  the  coordinate  system. 
The  analysis  is  predicated  on  knowing  the  initial  coordi¬ 
nates  of  points  A,  B  and  C.  Following  any  input  (9^),  the 
resulting  coordinates  of  A,B  and  C  are  calculated  via  a 
rotation  matrix. 

The  first  equations  of  constraint  are  those  of  constant 


link  lengths,  Figure  3.1. 

2  2  2  2 
(!„)  “  (Dx-Fx)  +(Dy-Fy)  +(Dz-Fz)  (3.1) 

UDE)2  -  (Dx-Ex)2+(Dy-Ey)2+(Dz-Ez)2  (3.2) 

UEF)2  -  (Ex-Fx)2+(Ey-Fy)2+(Ez-Fz)2  (3.3) 

(Aad)2  ■  (Ax- Dx) 2+( Ay- Dy ) 2+( Az-Dz) 2  (3.4) 

(*BE)2  "  (B*-Ex>2+(By-Ey)2+(Bz-Ez)2  (3.5) 

UCF)2  -  (Cx-Fx)2+(Cy-Fy)2+(Cz-Fz) 2  (3.6) 

UGD)2  -  (Gx-Dx)2+(Gy-Dy)2+(Gz-Dz)2  (3.7) 

(*GE)2  ■  (Gx-Ex)2+(Gy-Ey)2+(Gz-Ez)2  (3.8) 

(Arp)2  -  (Gx-Fx) 2+(Gy-Fy) 2+(Gz-Fz) 2 


(3.9) 


For  static  equilibrium,  Figure  3.2 


z  i?5  -  0 


£ADeA/D  +  fBE*B/E  +  £  CF*C/ F  "  “gj 


(3.10) 


E  M  -  0 
G 


IT,  x  f  1  +Tf.  xf  e  . 
D/G  AD  A/D  E/G  BE  B/E 


(3.11) 


+  \/G  X  £CF*C/F  "  ° 


When  the  unit  vectors  ere  introduced  Into  Equation  3.10  end 
tbe  multiplications  carried  out,  these  equations  become: 

£  f  £ 


-(Ax-Dx)f  *E( Bx-Ex )-t-r—  (Cx-Fx) 

xad  xbe  xcf 

£  f  £ 

-0 

:  i 

(3.12) 

~£(  Ay- Dy  )+7^i(  By-  Ey  )-H~(  Cy- Fy  ) 
*AD  *be  *cf 

f  f  f 

-mg  :  J 

(3.13) 

~^(Az-Dz)+TH(Bz-Ez)-l>T^(Cz-Fz) 
XAD  XBE  XCF 

-0 

:  k 

(3.14) 

Taking  the  cross  products  in  Equation  3. 

11  yields 

the  £ol- 

1  owing 
£ 

AD 


( Dy-Cy ) (Az- Dz)  -  (Dz-Gz)  (Ay-Dy)j 

+  -r^$(Ey-Gy)(Bz-Ez)  -  (  Ez- Gz  )  (  By- Ey  )2 
BE  t 

+  y^[(Fy-Gy)(Cz-Fz)  -  < Fz-Gz  )  (  Cy- Fy  )]  « 


0  :  l 


(3.15) 


x 

Y —  i(Dz-Gz)  (Ax-Dx)  -  (Dx-Gx)  (Az-Dz) ^ 


AD 


+  1(  E*-Gi5 )  (Bx-Ex)  -  (Ex-Gx)  (Bz-Ez)l 

*BE  * 

i. 


+  -~-|(Fz-Gz)  (Cx-Fx)  -  (Fx-Gx)(Cz-Fz)\ -0  :J 


(3.16) 


-^^(Dx-Gx)  (Ay-Dy)  -  (Dy-Gy)  ( Ax-Dx)^ 

*AD 

+  -j— — ■  i(ExrGx)(By-Gy)  -  ( Ey- Gy )  ( Bx- Ex  jl 
XBE 

t  _ 

+  —— A(Fx-Gx)  (Cy-Fy)  -  (Fy-Gy)  (Cx-Fx)l  -0  sic 


(3.17) 


Equations  3. 1-3. 9  and  3.12-3.17  are  then  solved  numeri¬ 
cally  for  the  coordinates  of  points  D ,  E  and  F  on  the  output 
link.  The  output  rotations  are  then  calculated  using  Rodri¬ 
gues'  formula  [2]  to  find  the  rotation  of  the  normal  to  the 
output  link  from  the  initial  assembly  position.  Rodrigues' 
formula  is  used  rather  than  a  general  displacement  matrix  to 
reduce  computational  time. 


3.3  Comparison  of  the  Spatial  and  Planar  Results 


Fint,  Rodrigues  formula  is  only  valid  when  a  body  is 
rotating  about  a  fixed  axis  without  any  trana la tl ons .  This 
was  verified  oy  calculating  the  displacements  of  the  output 
link  using  a  general  displacement  matrix,  and  then  noting 
that  the  displacement  terms  were  equal  to  zero. 


The  primary  results  of  this  analysis  are  displayed  in 

6 

Figure  3.3,  where  K'  is  plotted  against  —  for  different 

la 

values  of  The  ratio  K'  is  the  magnitude  of  the  rotation 

la 

of  the  output  link  about  the  X  and  Z  axes  over  the  magnitude 
of  the  input  rotation.  The  3  dimensional  analysis  of  the 

[ 2 ]  Bisshopp ,  K.E.,  "Finite  Rotations  of  a  Rigid  Body  and 
the  Airplane  Landing  Gear  Problem",  Jnl.  Mechanisms, 
Vol.  3,  Pergamon  Press,  1968  ,  pp.  203-208. 


2 


mechanism  indicates  a  smaller  K  than  that  predicted  by  the  2 
dimensional  analysis  due  to  rotations  on  the  order  of  hun¬ 


dreds  of  arc-seconds  of  the  output  link  about  the  Y  axis. 

e  h 

Figure  3.4  Illustrates  K'  verses  — ,  for  various  — 

ratios,  and  a  comparison  with  Figure  2.5.  As  previously 

shown  in  Figure  3.3,  the  spatial  analysis  yields  a  lower  K 

than  does  the  planar.  It  can  be  seen  that  the  difference 

between  the  spatial  and  planar,  at  each  j-,  is  similar  to 

that  seen  in  Figure  3.3,  and  that  the  curves  for  the  spatial 

h  e 

and  planar  models  begin  to  coalesce  as  ^  or  -  go  to  intin- 


Slnce  the  two  models  predict  the  same  performance  as  ~ 

la 

h 

goes  to  zero  or  ~  goes  to  infinity,  we  can  conclude  that  the 

la 

planar  analysis  is  a  suitable  substitute  for  Che  spatial 
model,  with  the  understanding  that  K  found  from  the  planar 
is  larger  than  K'  found  from  the  spatial  analysis. 


Figure  3.5  illustrates  the  change  in  one  of  the  output 
link  lengths,  while  the  two  others  are  held  equal.  The  fig¬ 
ure  shows  that  K'  actually  decreases  below  the  nominal 
value,  all  lengths  (L)  equal,  as  one  length  is  increased  to 
roughly  1%  larger  than  the  others.  Let  the  input  rotation 
be  about  the  axis,  Illustrated  in  the  inset  of  Figure  3.5, 


which  is  measured  from  the  positive  X  axis  parallel  to  the 
link  which  is  changing,  described  by  the  angle  a.  The  lower 
curve,  for  which  the  change  in  K'  is  very  severe,  is  the 
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function  of  n/L  for 
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CHAPTER  4 


PLANAR  MODEL  ACCOUNTING  FOR  HIRE  STIFFNESS 


4 . 1  Obj  actives 

As  previously  mentioned,  links  2  &  4  are  made  of  cable 
with  less  than  Ideal  joints  at  the  ends.  They  have  been 
assumed  to  have  Infinite  flexibility  so  that  the  bending  m 
the  cables  occurs  at  an  infinitesimal  area  at  the  the  con¬ 
nection  to  the  plates. 

For  this  section  the  cables  will  have  some  finite 
stiffness  and  are  subjected  to  shearing  forces  and  moments 
applied  to  the  ends.  A  planar  model  of  the  mechanism  is 
developed  here  which  Includes  the  stiffness  in  the  cables  as 
a  system  parameter.  It  will  be  shown  that  including  this 
effect  in  the  kinematic  representation  of  the  mechanism 
leads  to  results  which  will  agree  well  with  the  experimental 
data  gathered  in  the  next  section. 


Development  of  Eauetlona 


Tbe  followii 


lumptioas  bare  been  nede: 


*  negligible  streln  energy  is  stored  in  tbe  cables  et  tbe 
essembly  position. 


*  Tbe  cables  ere  Eu ler-fiernou 1 1 i  beans  sucb  tbet  tbe  cur* 
veture  is  proportional  to  tbe  bending  moment. 


*  Deflection  of  tbe  beams  is  smell  sucb  tbet  bending  does 
not  alter  the  length  of  tbe  cables. 


'  *  Strains  in  tbe  individual  wires  of  tbe  cable  ere  small 
and  tbe  axial  and  bending  responses  are  additive. 

*  Friction  and  twisting  in  tbe  cables  are  neglected. 

4.2.1  Force  and  position  solution 

Upon  assembly  of  tbe  mecbenlsm,  the  cables  are  plasti¬ 
cally  deformed  around  tbe  junction  at  tbe  input  link  and 
subsequently  assume  an  unstressed  kinked  shape  at  tbls  point 
if  disassembled.  Tbe  assumption  of  no  initial  stored  strain 
energy  is  then  an  approximation  to  tbe  state  of  the  cables 


at  assembly.  This 


iptlon  will  lead  to  a  pinned  boundary 


condition  at  tbe  input  link  whereas  a  clamped  boundary  con¬ 
dition  will  be  used  et  the  output  link.  Tbe  cables  will  be 
constrained  to  make  small  deflections  about  the  assembly 
position  in  order  to  simplify  tbe  analysis.  These  assump¬ 
tions  lead  to  a  schematic  representation  of  tbe  mecbanlsm  as 
illustrated  in  Figure  4.1. 


Figure  4.2 


Dlegrem  of  moments  end  forces  acting  on 
the  outout  link. 


L«t  6  tad  6  be  nttiurid  from  the  initial  conf izurt- 
A  B 


tion  of  the  ■icbtolaa  upon 


ily.  6  and  0  auat  ba 
A  B 


altered  by  the  rotation  of  the  input  link,  0^  ,  in  order  to 


maintain  the  axial  force  parallel  and  the  shear  force  per¬ 
pendicular  to  the  undeforsed  cable.  The  subsequent  angles 


after  a  rotation  of  0.  will  be  and  as  shown  in  Figure 

X  Ad 


4.2. 


*A  “  9A  '  9i 


♦«  “  9„  +  © 

B  B  1 


For  static  equilibrium,  the  force  and  moment  balance  on 
the  output  link  yields. 


E  F„  -  0 
3 


Ha«s^  +  PASin*A  +  HflCos^B  -  PBsin$B  -  0  (4.1) 


P  cost  -  H  sine  +  P  coit  +  H  slnt  -  mg  -  0  (4.2) 

AAAABddB 


Em  -  o 

A 

-“a  '  MB  +  PBLCOs(n"0o>  +  HBLsln(*B-0o) 

L 

-  mg  -  6)cos0q  +  esin0o  »  0  (4.3) 

Referring  to  Figure  4.3,  the  deflections  of  the  cables 
are  along  a  line  perpendicular  to  the  undeforsed  initial 
position,  so  that  point  A  can  be  calculated  from 


A 


W  cos* 


Similarly, 


Ay  -  -  V*A  +  WAginn 

Bx  "  6c°a0i  +  V10*B  ‘  WBC°*S 

By  -  6.in01  -  Abco.«b  - 

T fie  length  of  the  bottom  plate  remains  unchanged,  thereby 

connecting  the  deflections  w  and  w  . 

A  o 

(Bx-*x)a  +  <By-*y>2  -  t/  (*•*) 

B  -  A 

—2 - -X-  -  tan0  (4.5) 

B  -  A  o 

x  x 

4.2.2  Derivation  of  beam  equations 

The  cables  have  been  assumed  to  be  beams  which  are 
loaded  primarily  in  tension.  The  Important  equations  will 
be  derived  for  one  cable  with  the  second  cable  being  loaded 
in  a  similar  way. 

Referring  to  Figure  4.4,  and  taking  moments  about  point 
p: 

N  -  M  +  P(w.  -  w)  -  H,(Jt  -  u)  -  0 


Y 


Deformed- 

Cable 


Undeformed 

Cable 


Figure  4.3  Deflection  diagram  for  both  cablea  aa  uaed 
for  the  poaition  aolution  of  pointa  A  and  B 


Figure  4.4 


Free  body  diagraaa  for  one  cable  indicating 
the  variablea  in  the  beam  equationa. 


Sloe*  to*  cables  have  bten  assumed  Co  be  Euler  Bernoulli 


Deans  with  small  curvature,  we  can  then  write: 


2 

d  w 


M 

El 


<1  M  2  2  P 

- — —  •  P  N  -  0  where  P  ■  —  (4.6) 

.  2  r  u  El 

du 

where  El  Is  the  flexural  rigidity  of  the  cable.  The  general 
solution  of  Equation  4.6  is 


or 


-  C^coshpu  +  sinhpu 


- j  -  (C^coshpu  +  C^slnhPu)  (4.7) 

duZ 

As  previously  discussed  in  section  3.2.1,  the  boundary 
conditions  will  make  use  of  the  fact  that  the  cables  are 
nearly  pinned  at  the  Input  link  and  nearly  clamped  at  the 
output  link.  The  boundary  conditions  can  be  visualized  with 
the  aid  of  Figure  4.5  where  the  cable  and  plates  are  magni¬ 
fied  to  show  the  connections.  The  cable  Is  firmly  clamped 
at  the  bottom  plate  via  a  screw  and  washer.  The  cable  at 
the  top  is  drawn  around  a  knife  edge  yielding  a  joint  which 
has  the  aspect  of  being  sensitive  to  plate  rotation  while 
not  exerting  a  moment  on  the  plate.  So  a  pin  Joint  is  a 


reasonable  approximation  to  the  actual  boundary  conditions 


The  cables  are  not  securely  fastened  to 
the  input  link  at  the  knife  edge. 

Since  the  attachment  is  away  from  the 
knife  edge,  a  pinned  condition  occurs. 


The  cables  are  securely  clamped  to 
the  output  link  by  screws. 


Figure  4.5 


Sketch  of  the  Instrument  illustrating 
the  cable  boundary  conditions. 


at  the  top  link 


Tb«  boundary  conditions  can  then  be  stated 


as  : 


u  -  0 

u  -  i 


w  ■  0 

d  v 

-  © 

d  u  o 


2 

d  w 
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d  w 


The  above  boundary  conditions  lead  to  the  evaluation  ot 
tbe  constants  in  Equation  4.7. 


C 

1 


c 

4 


0 


M 

c _ i_ 

2  sinhfH 


C3  "  0o  + 

Then  the  displacement  and  shear  at  the  bottom  of  the  cable 
can  be  found  from  equation  4.7  and  its  derivatives. 

■  -^(C^cosupA  +  C2sinnpjO  +  C^X  +  C4 
Inserting  the  constants  yields 

•*  '  «.  +  -  »  <*••> 
From  the  equation  for  the  shear  on  the  cable 


ftMcoshpf 

Pslnhpf 


El 


— (  C  ^  sinti|3Jl  +  C^coshpA) 


ir 


-jj—  -  C^slnnpi  +  C^coshfJX 


Inserting  the  constant!) 


H  .  - 


COShpA 

X  s lnh£A 


(4.9) 


Equations  4.8  and  4.9  are  then  written  for  each  cable. 
With  the  Inclusion  of  equations  4.1  -  4.5  there  are  enough 
constraint  equations  to  solve  for  the  equilibrium  position 
of  the  output  link  numerically. 


4.3  Results  from  stiffness  analysis 

In  order  to  evaluate  the  magnitude  of  the  stiffness 

effects  on  K,  the  moment  of  inertia  (I)  and  the  modulus  of 

elasticity  (E)  were  found  for  a  0.006  in.  diameter  7X1 

stainless  steel  cable.  The  moment  of  inertia  for  the  cable 

- 11  4 

was  calculated  to  be  4.494  x  10  In  ,  while  the  modulus  of 

6 

elasticity  was  determined  to  be  approximately  14  x  10  psl 

6 

[ 3 ] [ 4 ]  .  Stainless  steel  has  a  modulus  of  28-29  x  10  psl, 
but  a  correction  has  been  made  for  slippage  and  the  reduced 
area  of  the  cable  as  opposed  to  a  solid  wire.  These  values 
are  only  approximate  as  no  experimental  results  were 

[ 3 ]  Costello,  G . A . , e tal , "Ef f ec ti ve  Modulus  of  Twisted  Wire 
Cables",  Journal  of  the  Engineering  Mechanics  Division, 
ASCE,  Vol.  102,  No.  EMI,  February,  1976,  pp.  171-181. 

[4]  Shlgley,  J.E.  ," Mechanical  Englneerlug  Design",  Tnird 
Edition.  McGraw-Hill,  New  York,  1977. 


obtained  to  confirm  the  moment  of  inertia  or  modulus  of 


elasticity  of  the  cable. 


Figure  4.6  illustrates  the  influence  of  the  stiffness 

1  6  h  h 

on  ~  as  a  function  of  —  and  ~.  The  effect  of  —  is  shown  by 

successive  plots  starting  with  -  1/2  (Figure  4.6a)  and 

increasing  up  to  —  ■  10,  (Figure  4.6d).  As  with  the  earlier 

la 

h 

analysis,  ~  >  1  has  little  effect  on  the  input  output  ratio 
K.  In  fact,  there  is  no  difference  between  Figures  4.6c  & 
d  . 


6 


The  ~  ratio  is  plotted  on  each  figure  tor  what  klnemat- 
L 

6 

ically  would  be  a  pure  pendulum,  rj-  -  0 .  K  does  not  go  to 
zero  as  it  would  in  the  kinematic  analyses,  but  Instead  is  a 

strong  function  of  El.  It  is  apparent  from  Figure  4.6  that 

6 

the  kinematic  Influences  at  small  —  ratios  are  negligible 

6 

compared  to  the  stiffness  Influences.  As  —  increases  the 

kinematic  errors  begin  to  out  weigh  the  stiffness  errors. 

Therefore  the  earlier  analyses  are  valid  for  small  values  of 

6 

El  or  for  large  values  of  — . 


e 

In  Figure  4.7  the  effects  of  —  are  illustrated.  Fig- 

L 

,  h  6 

ures  4.7a  through  4.7d  are  based  on  ~  ■  1  while  —  is  lncre- 

e 

mented  on  each  plot.  Again  —  has  a  strong  influence  on  K 

and  and  to  minimize  the  positioning  error  of  the  mechanism  a 

,  e 

large  —  is  recommended. 


In  Chapter  2  it  was  shown  kinematically  that  the  CG 
could  be  located  anywhere  up  to  50%  of  L  along  the  output 


link  without  an  associated  error.  The  Inclusion  of  stiff¬ 
ness  into  the  kinematic  model  shows  that  A  may  vary  widely 
without  introducing  any  errors  and  also  that  the  maximum 
limits  are  dependent  on  the  cable  stiffness. 


A 

Figures  4.8a, b  illustrates  of  30%  and  20%  reapec- 

A 

tively.  The  curves  do  not  differ  from  that  of  —  ■  0,  indi¬ 
cating  that  mlscentering  of  the  CG  is  not  as  critical  as 
indicated  in  Chapter  2.  Because  the  bending  of  the  cables 
is  creating  a  reactive  moment  on  the  output  link  the  sign  of 
A  is  important.  If  A  is  Increased  in  the  direction  of  the 
deflection,  (l.e.  to  the  left  in  Figure  4.4)  the  maximum 

value  for  £  is  approximately  30%.  On  the  other  hand,  if  A  is 
negative  with  respect  to  the  displacements,  the  maximum  is 
close  to  40%  before  the  cables  are  loaded  in  compression. 

These  results  are  dependent  on  El  and  are  stated  for 
2 


-  1000, 


or  the  first  point  plotted  in  Figure  4.8. 
these  values  would  increase  to  50%. 


As 


CHAPTER  5 


EXPERIMENT 


5.1  Op j  ec  tives 

Experimental  results  were  needed  to  Identify  problems 
with  the  mec&anlsm  end  to  give  e  bench  mark  tor  the  analyt¬ 
ics.  Three  mechanisms  were  tested,  two  with  a  large  Input 
range  and  one  for  much  smaller  Inputs.  It  was  desired  that 
these  tests  Indicate  the  roles  of  the  geometry  of  the 
mechanism,  hence  three  different  geometries. 

5.2  Description  of  Experiments 

The  experimental  procedure  was  divided  Into  a  small 
Input  teat  (for  Input  motions  from  10  to  300  arc-seconds) 
and  a  large  Inpuf  test  (for  rotations  over  300  arc-seconds). 

The  small  Input  test  is  shown  In  Figure  5.1a.  An  auto- 

colllmator  was  used  to  monitor  both  input  (@t)  and  output 

(9  )  by  use  of  a  beam  splitter.  The  autocollimator  enabled 
o 

a  0.5  arc-second  accuracy  to  be  held  on  both  the  Input  rota- 


Figure  5.1  Sketch  of  epperetue  used  to 
determine  K  experimentally. 
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Since  cite  au  toco  1  lima  to  r  bad  only  a  ±  5  arc-minute 


range,  the  Input  rotations  during  the  coarse  experiment  were 
monitored  by  using  a  spirit  level  wit n  5  arc-minute  gradua¬ 
tions  as  shown  in  Figure  5.1b  Experimental  procedures  con¬ 
sisted  of  rotating  the  lab  stand  to  a  predetermined  0  via  a 

it 

micrometer  and  then  reading  the  output  angles  by  use  of  the 
au  toco  1 lima  tor . 

5.3  Results  from  Static  Positioning  Experiments 

Tests  were  conducted  on  three  mechanisms  which  are 
described  by  the  table  (Table  5.1)  of  non- d imens ionless  link 
lengths.  The  major  differences  between  the  mechanisms  were 
the  parameters  (6/L).  Both  of  the  cases  A  and  B  have  rela¬ 
tively  small  input  links  (6)  compared  to  case  C.  The  height 
(h)  of  A  and  B  was  again  much  smaller  than  C  as  indicated  by 
the  h/L  ratio  in  the  second  column.  Plus,  the  stiffness 
mgf2 

parameter  - g-- —  for  case  A  was  much  smaller  than-  cases  B  and 
t  I 

C  due  primarily  to  an  increased  mass  of  these  mechanisms. 


Table  5.1 

Non-Dimensional  Ratios 
Describing  Experimental  Mechanisms 


n 

6 

L 

h 

L 

A 

L 

e 

L 

E9 

A 

0.0034 

0.722 

0.0 

0. 144 

9342 

B 

0.0021 

0.650 

0.0 

0.072 

12160 

C 

0.0197 

1.559 

0.0 

0 . 144 

12160 

5.3.1  Small  Input  test 


Eight  experimental  states  are  plotted  In  Figure  5.2  for 
the  Case  A  mechanism  along  with  the  theoretical  values  from 
the  planar  kinematic  and  planar  stiffness  models.  At  each 
Input  state,  an  error  band  Indicates  the  range  of  the  meas¬ 
urements  while  the  circle  Indicates  the  average.  Linear 
regression  yields  the  best  fit  line  through  the  data  and  the 

slope  of  that  line  Is  K  .  Comparisons  of  K  to  the 

exp  exp 

kinematic  model  of  Chapter  1,  K  and  the  kinematic  model 

kin 

Including  cable  stiffnesses,  K  are  then  shown  In  Table 

stlf 

5.2. 

It  is  apparent  from  Table  5.2  that  the  kinematic 

analysis  yields  a  K  which  is  extremely  low  In  comparison  to 

the  experimental,  with  a  difference  of  78. 6%  from  K,  .  to 
r  kin 

K  .  The  second  analysis  gives  a  much  better  lndlcetor  of 
exp 

the  performance  of  the  mechanism  with  only  a  13.4%  differ¬ 
ence  between  K  and  K  , .  The  primary  reason  for  the 

exp  stlf  7 

large  change  between  K.  ,  and  K  is  due  to  the  small 

kin  stlf 

value  of  As  seen  in  Chapter  4  Figure  4.6,  the  effects  of 

L 

6 

the  stiffness  become  predominant  as  —  becomes  small. 


K 

%  Diff 

s  tlf 

0.0162 

13.4% 

0.0144 

-16.1% 

0.0221 

-13.9% 

%  Dif  f 


0.0187  0.0040 

0.0124  0.0020 


0.0194 


0.0185 


K 


(%  Diff.  defined  as 


5.6% 


analysis 


X  100%  ) 


5.3.2  Large  input  test 


The  results  from  the  large  Input  test  for  case  B  are 

Illustrated  by  the  lower  set  of  curves  in  Figure  5.3.  Case 

B  consistently  has  a  smaller  angular  error  ( 0Q )  than  does 

Case  C  as  predicted  by  the  planar  kinematic  analysis  of 

each.  The  kinematic  analysis  again  predicts  a  lower  K  than 

that  found  experimentally,  as  Indicated  in  Table  5.2  with  a 

85.6%  difference  of  K  to  K  .  Although  the  stiffness 

kin  exp 

model  reflects  the  experimental  data  better,  with  a  16.1% 

difference  from  K  to  K  ,  it  also  predicts  a  larger 

stiff  exp 

value  for  K  than  that  measured.  This  is  surprising  in  view 

of  the  case  A  results  where  K  was  larger  than  K 

kin  exp 

Turning  to  case  C,  Illustrated  by  the  upper  set  of 
curves  In  Figure  5.3,  It  is  apparent  that  the  planar 
kinematic  analysis  for  this  mechanism  accurately  predicts 
the  experimental  results  with  a  difference  of  only  5.6%  as 


test  for 


shown  In  Table  5.2 


This  shows  that  the  kinematic  errors 


due  to  the  large  —  are  outweighing  the  errors  due  to  stiff- 

la 

ness.  The  Inclusion  of  stiffness  in  the  analysis  Increases 

K  ....  to  0.0221,  for  a  difference  of  13.9%  between  K,  and 
stiff  kin 

K  ^  with  the  stiffness  model  again  yielding  the  larger  K. 

5.3.3  Discussion  o f  results  and  error  sources 

As  mentioned  previously,  the  symbols  on  the  experimen¬ 
tal  curves  reflect  the  average  of  the  data  at  that  point, 
while  the  error  bands  extend  above  and  below  these.  The 
error  band  was  as  much  as  ±  30%,  15%,  and  10%  of  the  average 
values  for  cases  A,  B  and  C  respectively.  For  case  A  the 
large  spread  is  partially  due  to  the  resolution  of  the  auto¬ 
collimator  (i  .5  arc  seconds)  while  a  1  arc  second  measure¬ 
ment  was  being  made.  Hence,  an  accuracy  problem  exists. 
For  cases  B  and  C  the  spirit  level  used  to  monitor  the  input 
could  be  read  to  ±  1  arc  minutes,  while  the  magnitude  of  the 
input  angles  were  on  the  order  of  1  degree.  So  it  can  be 
concluded  that  cases  B  and  C  were  run  with  reasonable  accu¬ 
racy,  while  the  data  for  case  A  could  be  different  from  the 
nominal  values  plotted  in  Figure  5.2  or  listed  in  Table  5.2. 

Since  the  stiffness  analysis  predicts  K  values  greater 
than  the  experimental  for  cases  B  and  C,  it  indicates  that 
the  analytics  are  incorrect  or  that  some  variable  or  values 
were  incorrect.  As  mentioned  in  Section  3.3,  the  value  of 
El  was  merely  an  estimate.  If  El  were  decreased  by  a  small 


V  *»  *•  \  ",  •„  *.  !*.*»*.' 

»  ,*  ,*  ,*  ( 


54 


mg  Jl  3 

amount  K  would  decrease  sizably  since  — gy-  ■  14  X  10  is 
seen  to  lie  on  the  rapidly  increasing  portion  of  the  curve 

shown  in  Figure  3.6b.  From  this  figure  one  sees  that  a 

6 

decrease  of  E  from  14  to  12  X  10  ,  a  decrease  of  14%,  would 
change  K  approximately  6%.  Also,  the  difference  between  the 
planar  and  spatial  analyses  from  Chapter  3  was  roughly  10%, 
thereby  reconciling  the  differences  between  the  experimental 
results  and  the  results  from  the  kinematic  model  including 
stiffness  as  being  dependent  on  the  value  tor  E  and  the 
rotations  out  of  the  plane  of  the  planar  analysis. 

It  should  be  noted  that  errors  could  be  introduced  in 
the  planar  analysis  including  stiffness  where  small  deflec¬ 
tions  were  assumed.  Inherent  in  this  assumption  is  that  the 
distance  along  the  undeformed  cable,  (u)  from  Chapter  4  Fig¬ 
ure  4.7,  remains  constant  and  the  deflections  (w)  are  from 
the  undeformed  cable.  The  actual  length  of  the  caole  is 
1~ 2  2 

A  »  \|u  +  w  .  While  this  is  a  small  effect,  it  would 

change  K  by  1%  for  example  in  Case  C.  So  a  more  accurate 

analysis  would  be  to  assume  large  deflections  of  the  cables. 

The  boundary  conditions  assumed  for  the  cables,  how¬ 
ever,  are  the  most  representative  of  tne  actual  situation. 
If  a  clamped  condition  were  to  be  assumed  at  the  input  link, 
the  system  would  be  stlffer  and  raise  K,  while  a  pinned  con¬ 
nection  at  the  bottom  would  lower  K,  yet  be  unrealistic  due 

to  the  firm  attachments.  Therefore,  in  comparison  with  the 


CHAPTER  6 


DAMPER 


6 . 1  Objectives 


As  previously  mentioned,  the  Level  will  be  used  in  an 
environment  in  which  vibrations  are  present  due  to  the 
machine  tool  the  instrument  is  monitoring  and  others  in  the 
surrounding  production  area.  These  vibrations  must  be  elim* 
lnated  by  a  damper  because  the  leveling  mechanism  is  passive 
and  has  no  means  o£  compensation.  The  characteristics 
looked  for  in  the  damper  were:  1)  the  ability  to  function 
equally  well  in  two  orthogonal  directions,  2)  remove  oscil* 
lations  above  1%  of  the  input  (100  arc  seconds  nominally) 
and  3)  not  interfere  with  the  function  of  the  level. 


These  constraints  lead  to  tne  cylindrical  damper  shown 
in  Chapter  1  Figure  1.4,  and  again  in  Figure  6.1.  Air  will 
be  used  as  the  viscous  medium  not  because  it  is  the  optimum, 
but  because  seals  are  not  required  and  nence  no  friction 
will  be  Introduced  into  the  system  to  degrade  the  mechanisms 
positioning  accuracy.  The  analytics  developed  herein  will 
be  used  to  determine  the  feasibility  and  then  the  design  of 
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a  cylindrical  air  damper  for  this  application. 

6.2  Cylindrical  Pandulum  Daaper  Mecbanlcs 

Equation  6.1  [5],  the  Reynolds  Equation,  defines  the 
mecbanlcs  of  a  fluid  between  walls  Inclined  to  each  other 
and  moving  relative  to  each  other.  It  is  assumed  that  the 
space  between  the  walls  is  small  enough  to  suggest  that  the 
flow  is  laminar  and  fluid  velocity  perpendicular  to  the 
walls  is  negligible  compared  to  that  parallel  to  the  walls. 
These  conditions  prevail  approximately  if  the  ratio  between 
the  developed  flow  path  length  and  the  gap  is  better  than 
1000:1.  The  required  ratio  is  even  larger  if  the  fluid 
viscosity  is  very  low  as  is  the  case  with  air.  As  a  begin¬ 
ning,  approximations  based  on  laminar  flow  conditions  will 
be  assumed  even  though  the  length-gap  ratio  may  not  be  quite 
optimum. 

Reynolds  Equation: 


<VA 

•X"!' 


f(x,y)  -  -  U- 
See  list  of  symbols  for  Chapter  6. 
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5y 

’  e  *  - 

5  h 

dh 

(6.1) 

,  •' « ■ 

'ftx 

W-r-  +  2  V 
dy 

W- 

The  Reynolds  Equation  can  be  written  in  polar  coordi¬ 
nates  consistent  with  the  damper  geometry.  Further  assuming 
constant  viscosity  of  the  air 

( 5  j  Cameron ,  A. /'Principles  of  Lubrication",  J.  Wiley  Publ. 


at  any  level  y  the  gap  h  can  be  written  as 


h  -  c  +  • ( y ) cos©  or  a  -  c  -e(y>cosP 
To  get  a  general  estimate  of  the  mechanism  behavior,  it  Is 
advisable  to  non-dimenslonallze  Important  parameters.  For 
the  pendulum  damper,  the  following  transformations  are  con¬ 
venient. 


H(y) 


h(y ) 
h 
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e(y) 


«(y) 


h  (y-A)  -  h  (y-0) 
h  (y-0) 


[H  -  »1 
l-T - 2 1 

l_  o  J  p  region 


h  (y-o)  -  h  (y-A) 

b  - - iTT^-TT - 


region 


In  a  pendulum  damper,  the  U  and  W  velocities  are  inef¬ 
fectual  compared  to  the  squeeze  component  V  so  that  Equation 
6.2  can  be  simplified  and  with  the  non-dimensional  transfor¬ 
mation  becomes 
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vjr 
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-0  +2V] 


The  above  partial  differential  equation  can  not  be 
solved  in  a  closed  form  but  a  simplification  suggests  Itself 


T 


by  the  reelitles  of  tilting  dampers.  Because  h  varies  along 

y  it  is  not  practical  to  use  long  damper  cylinders.  If,  on 

,  *  1  . 

the  other  hand,  short  cylinders  are  used  (  “  <  ~  )  it  might 
be  assumed  that  pressure  gradients  in  the  circumferential 
direction  will  be  considerably  less  than  in  the  axial  direc¬ 
tion.  If  that  is  the  case,  the  first  expression  on  the  left 
hand  side  of  the  Reynolds  Equation  might  be  ignored  reducing 
the  partial  differential  equation  to  the  ordinary  differen¬ 
tial  equation  6.3. 


where 

as 


dY  3 

n 

o 

the  squeeze  velocity  around  the  cylinder  is 


(6.3) 

described 


$(L-y)cos9 


$L(  l-Y-jr)cos0 


and 


V  -  <frL(  1-Y-£)cosP 


therefore 


dY 


l-Y-|)cosp 


Let 
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»o  that 


/  3dP, 
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dY 
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Since  Y^  <  1  let  ua  expedite  the  search  for  behavior 
characteristics  by  setting  the  right  hand  side  to 

X 

(1  -  0)cos(3  .  A  refined  analysis  with  Y~  included  might  be 
pursued  later. 


Reduced  Reynolds  Equation: 


,  3dPv 
d ( H  ) 


dY 


-  COS  P 


(6.4) 


At  any  position  p  integrating  6.3  yields 
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dY 
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Integrating  again 
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(6.5) 


From  the  non-d imens Iona  1  transformations  in  the  p  region 


H  -  1  +  aY 


therefore 
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dY 
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Substituting  the  above  into  Equation  6.4  yields 
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or 
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AC  may  position  p  a  strip  RdP  wide,  of 
cylinder,  has  a  force  dP  exerted  on  It: 

dF  ■  / pRdpdy  ■  /PdY 

0  °  0 


dF 


H 

Rlcospdp  j 
a3D  H 

o 


H  H 

I+  (1  +-^)xo 


RlcosPdp  f  ,  /  .  .  v  .  a  "1 
it - — 

a  D  *■  — 

The  force  component  In  the  direction  of  the 

change  Is 


Ricos 


JMf. 


(6  .6) 

the  damper 
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eccentrlcl ty 
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1.(1  +  .)  +  7-i— 1 


In  the  region  -w  <  p  <  —  the  total  force  generated  la 


« 

T  2 


where 
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0  a  u  J 
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“7^  ‘ 


AH(p)  ■  Jt$cosp  &  HQ(p)  ■  1  -  e^cosp 
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Therefore 
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A  similar  treatment  to  the  9  region  yields  an  analogous 
result  with  suction,  rather  than  positive  pressure,  being 
generated,  hence,  the  two  regions  reinforce  each  other  so 
that  the  total  fluid  force  becomes 


where 


F  -  24^VoR(-)  Kd 


(Te^)3  l  7^7e['(ln(l+*)  +  ln(l+b))  +  ITT  +  TTFj*0 


It  C  O  8  0 
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l( l-eQcos0) 
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It  CO80 
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L ( l+e^cos© ) 


ls  given  on  Figure  6.2  es  a  function  of  the  eccen¬ 
tricity  ratio  eo  and  length  ratioa  .  U  «  0  the  pendu¬ 
lum  is  exactly  centered  on  it's  cylindrical  housing.  At 
e  ”  1.0  the  pendulum  ls  touching  the  wall.  The  minimum 
damping  force  occurs  at  e  ■  0 


There  were  several  approximations  made  to  expedite  the 
analysis.  While  the  results  are  not  exact,  the  errors  are 
probably  sufficient  for  design  purposes  and  they  do  show  the 
role  of  the  Important  parameters. 

6.3  Pendulum  with  Fluid  Damping 

As  a  convenient  but  reasonable  approximation  to  the 
actual  3  dimensional  dynamics,  a  2  dimensional  analysis  was 
developed.  A  sketch  of  the  damper  system  is  shown  in  Figure 
6.3. 


Taking  moments  about  the  pivot  0  (a  fixed  Instantaneous 
center  of  curvature  rather  than  the  perturbated  small  move¬ 
ment  was  assumed): 


Hence 


£  -  mgLSint  -  CDLD$  “  J9 


J$  +  ♦  +  ®gLSin^  ■  0 

or  for  small  motion. 


♦  + 


♦  +  5*  L* 
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The  undamped  natural  frequency  of  the  pendulum  is: 


L 

:J 


The  damping  factor: 


L(l+(lc/L)2) 


C  L 
DP 

2  Joi 


C  to  L 
DnD 


c  - 


C  L 
D  D 


2m(L)1,5  \ | g ( l+( k/ L ) 2 ) 


The  damping  force  is  described  as  F^ 

Section  6.1  F  is  also  described  as 
D 


C^V  but 


H  ence 


Fd  -  24nRV(-i)3K0 


CD  -  24^(-i)3KD 


Substitution  of  this  into  Equation  6.9  yields 


c 


24uR( 1/c)3 


u>  L_ 
n  D 


*  Kn 


(6.7) 


(6.8) 


(6.9) 


from 


(6.10) 


(6.11) 


(6.12) 


6.4  Experiment 


Equation  6. 1.2  demonstrates  the  design  parameters  and 
their  role  for  acquiring  a  working  damping  coefficient. 
This  was  used  to  construct  the  experimental  apparatus  of 
Figure  6.4.  The  experiment  was  fairly  simple  in  that  the 
only  instrumentation  required  was  a  laser.  The  laser  was 
directed  to  a  mirror  attached  on  the  mechanism  and  then  on 
to  a  sheet  of  graph  paper  positioned  at  the  correct  height. 
By  virtue  of  the  distance  from  the  mirror  to  the  graph 
paper,  a  magnification  of  the  mirror  rotations  occur,  such 
that  up  the  laser  beam  was  displaced  about  6  cm.  for  a  one 
degree  swing  of  the  mirror.  Then  the  cusps  of  the  laser 
beam  were  marked  on  the  paper  as  the  time  was  recorded. 


The  experimental  damper  was  sized  as  shown  in  Figure 
6.5.  The  damping  ratio  was  then  determined  as  follows.  The 
clearance  (c)  is 

D-D 
o  1 

c  -  - - -  -  .0135 

Substituting  the  data  in  Figure  6.3  into  Equation  6.12 
yields 


C  *  24  *  0.25  x  10 
or 


x  1.87 5 ( 


.5 


.0135 


>3  x 


11.4 


2  x  0.72 


2.56 

2.66 


C  -  0.11  kd 


At  £  «  0,  K. 
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?  L 


(Figure6.2);  C 


014 


,C|i? 


two 


and  at  e  ■  0.7;  X  ”  ~ — rr  “  .2;  ■  0.85 

4  .30  U 

If  X  Is  Increased  to  1“  and  every  otter  property  is 
unchanged  then 

,  1.3  ,2.84.2 

5  “  (. 5 *  X  *2.50*  X  *LLKD 

and 

at  e  -  0  £  -  .14  at  e  -  0.7  £  -  .93 

6 .5  Experimental  Results 

The  first  experiment  was  conducted  with  X  “  0.5".  The 
response,  Figure  6.4a,  Indicates  £  ”  0.17.  A  second  experi¬ 
ment  with  X  m  1.0  yielded  the  response  of  Figure  6.4b  and 
£  ■  0.73.  The  experimental  set  up  did  not  allow  a  determi¬ 
nation  of  the  eccentricity  In  the  dormant  state,  further¬ 
more,  the  resolution  possible  was  coarser  than  e  *  0.3. 

Experimental  C  sad  theoretical  £  are  0.17  and  0.094 
respectively,  from  the  first  test  (A  ”  0.5").  Experimental 
and  theoretical  £  are  0.55  and  0.93,  respectively  from  the 
second  test  (X  -  1.0").  Theoretical  values  were  besed  on  a 
0.3  eccentricity  between  damper  sleeves  considering  that  the 
test  set  up  static  e  is  only  an  estimate  and  the  fact  that 
the  rapid  decay  of  oscillations  makes  it  difficult  to  meas¬ 
ure  the  second  and  third  decayed  amplitudes;  agreement  witn 
theory  seems  adequate  for  design  purposes. 


*  .  »  V  V-  ■_>  "4  *  *  ■  •/  •>•*•*•*  V* 


•*-  A  •*. 


/  <*  • 


la  possible  co  design  an  air  damper  which  will  have  a  set* 
Cling  cime  of  l  second.  Also,  the  size  of  Che  damper 
required  for  this  example  Is  reasonable.  The  size  of  Che 
damper  Is  not  prohibitive  given  Che  size  of  Che  output  link, 
both  nearly  3.75  inches. 


Equations  6.12  witb  Figure  6.2  should  serve  as  adequate 


design  guides.  From  Equation  6.11  it  is  apparent  that  to 
maximise  damping  the  following  parameters  should  be  con¬ 
sidered  . 

I 

1.  The  wetted  length  to  radial  clearance  ratio  —  is  very 
critical  coming  in  as  a  cubic  influence. 

2 

L 

2.  The  ratio  - —  calls  for  setting  LQ  as  long  as  possible 

L1*5 

and  U  as  short  as  possible.  By  the  kinematics  of  the 
suspension  it  is  suggested  that  maximizing  L  may  serve 
both. 

3.  Both  the  mass  and  to  a  lesser  degree  the  radius  of 
gyration  should  also  be  minimized. 


Finally,  one  might  ask  how  small  a  value  of  C  is  toler¬ 
able.  In  other  words,  how  long  will  It  take  the  pendulum  to 
subside  to  a  motion  1%  of  the  initial  disturbance.  This 
decay  rate  is  approximated  by 


n 


1 

Tin  100 


In  100 

2*C 


n  -  number  of  oscillations. 


If  the  period  of  the  pendulum  is  .6  secs  and  £  »  . 1  the 
time  interval  is  4.4  seconds.  If  a  1  second  interval  is 
desirable  the  required  damping  factor  must  be  .44.  Hence  it 


Is  possible  to  design  an  air  damper  which  will  have  a  set¬ 


tling  time  of  1  second.  Also,  the  size 
required  for  this  example  is  reasonable, 
damper  is  not  prohibitive  given  the  size  of 
both  nearly  3.75  inches. 


of  the  damper 
The  size  of  the 
the  output  link. 


CHAPTER  7 


CONCLUSIONS 


A  proposed 

mechani 

sm  to 

level 

0 

p  tic 

s  a  bo 

u  t  two 

ho 

rizon- 

tal 

axes 

to 

a  high  d 

egree 

of  ac 

cu 

ra 

cy 

has 

been  s 

tud 

led  to 

determine 

the 

s  true  tu 

ra.l 

proper 

tl 

es 

that 

inf  lu 

enc 

e  the 

ope 

ration. 

Th 

e  mechan 

ism  1 

s  comp 

os 

ed 

o 

f  two 

terna 

*y 

links , 

one 

being 

the 

input  w 

1th 

the  o 

th 

er 

9 

the 

outpu 

t, 

be  ing 

s  us 

pended 

by 

mlnia  tur 

e  cab 

les  . 

Th 

e 

me 

chani 

sm  was 

de 

signed 

to 

simulate  a  a 

pherical 

pend 

ulum  s 

o 

tha 

t  ou 

tput 

ro  ta  tlons 

wou 

Id  be 

less 

than  1 

arc 

second 

f 

or 

i 

nput 

rotatl 

ons 

up  to 

100 

arc  se 

conds 

• 

Initially,  the  3  dimensional  mechanism  was  generalized 
to  the  two  dimensional  case.  A  kinematic  analysis  was  per¬ 
formed  the  two  dimensional  model  to  highlight  the  effects  of 
link  lengths  and  tolerances  on  the  positioning  errors.  From 
this  analysis  it  was  learned  that  the  ratio  of  the  output 
link  length  to  the  input  link  length  was  an  important  param¬ 
eter.  This  ratio  must  be  maximized,  up  to  1000:1,  to  get  a 
desirable  performance  from  the  mechanism.  The  vertical 
position  of  the  center  of  mass  was  also  an  important  system 
parameter  and  should  lie  below  the  bottom  link  to  ensure 


minimum  kinematic  error.  A  surprising  result  from  this 
analysis  was  that  the  center  of  mass  need  not  be  centered  on 
the  output  link  due  to  the  assembly  procedure  of  the  mechan¬ 
ism. 

Upon  comparison  of  the  results  from  this  analysis  with 
results  gathered  in  experiments,  it  was  found  that  the 
planar  model  predicted  lower  output  angles  for  two  of  the 
three  experimental  mechanisms.  It  was  hypothesized  that  the 
difference  was  due  to  rotations  of  the  output  link  which 
were  unaccounted  for  in  the  planar  model. 

A  spatial  kinematic  analysis  was  then  performed  which 
Included  all  of  the  geometry  found  in  the  experimental 
mechanism.  While  the  analysis  was  far  longer  than  the 
planar,  a  small  amount  of  additional  Information  was  gath¬ 
ered.  The  important  result  was  that  the  actual  3  dimen¬ 
sional  mechanism  has  less  positioning  errors  than  indicated 
by  the  planar  analysis  due  to  rotations  of  the  output  link 
about  the  vertical  axis  having  a  corrective  action  on  the 
rotations  about  the  horizontal  axes.  The  spatial  analysis 
checked  well  enough  with  the  planar  to  Indicate  the  former 
could  be  done  with  little  loss  in  accuracy.  Since  the  spa¬ 
tial  analysis  Indicates  lower  output  angles  than  the  previ¬ 
ous  analysis,  there  is  a  greater  discrepancy  with  the  exper¬ 
imental  data. 

The  effects  of  stiffness  in  the  mechanism  due  to  the 


beading  effects  of  the  support  cables  were  then  investigated 
to  understand  the  difference  observed  between  the  analytical 
and  experimental  results.  The  support  cables  were  thought 
to  be  slender  beams  pinned  at  the  input  link  and  clamped  at 
the  output  link  on  the  two  dimensional  representation  of  the 
mechanism.  This  analytical  approach  proved  to  accurately 
predict  the  output  observed  in  the  experimental  work.  It 
was  found  that  the  effects  of  the  moment  of  inertia  and 
modulus  of  elasticity  of  the  cables  were  stronger  parameters 
to  the  system,  for  certain  geometries,  than  the  kinematic 
properties.  It  was  shown  that  these  variables  increased  the 
output  angle  exponentially  such  that  the  minimum  modulus  and 
inertia  must  be  used  to  optimize  the  performance  of  the 
mechanism. 

An  air  damper  was  also  designed  to  minimize  distur¬ 
bances  of  the  optics  during  operation.  Two  cylinders,  one 
fixed  to  the  output  link  and  one  to  ground,  formed  a  two 
dimensional  damper  which  was  free  of  friction  while  dissi¬ 
pate  the  input  oscillations  effectively.  Again,  in  order  to 
establish  the  correct  geometry,  an  analysis  based  on  the 
Reynolds  Equation  assuming  constant  air  viscosity  was  per¬ 
formed.  Using  experimental  results  from  two  damper 
geometries,  it  was  shown  that  the  analysis  was  adequate  to 
characterize  the  system,  and  that  the  cylinder  length  should 
be  maximized  while  the  distance  from  the  input  link  to  the 
center  of  gravity  of  the  optics  should  also  be  maximized. 


RECOMMENDATIONS 


To  optimise  the  performance  of  a  two-axes  level  the 
following  should  be  considered. 

2 

mg 

1.  The  stiffness  ratio  should  be  maximized.  The  two 

possible  ways  to  do  this  are  reduce  the  cable  diameter, 
reduce  the  moment  of  inertia,  or  increase  the  wire 
length.  By  reducing  the  cable  diameter,  6  may  also  be 
reduced  and  by  increasing  the  wire  length,  the  height 
(h)  is  increased  both  yielding  better  performance  by 
reducing  the  stiffness  effects  and  kinematic  effects. 
The  mass  on  the  output  link  should  not  be  raised 
because  this  would  reduce  the  damping  ratio. 

2.  Maximize  the  height  (h)  of  the  mechanism.  As  illus¬ 
trated  in  Chapters  2  and  4,  the  h/L  ratio  is  the  next 
important  variable  after  6/L  and  the  stiffness  ratio. 
Certainly  one  would  not  elect  to  make  h/L  much  less 
than  unity,  as  K  increases  rapidly  after  this  point. 
The  damping  ratio  is  lnversly  proportional  to  the 
height,  Equation  6.12,  while  K  runs  ex poden tia 1 ly  with 
h.  Figures  4. 1,4. 2  and  4.3,  so  one  must  comprlmise  the 
settling  time  for  the  positioning  error,  but  the 
minimum  h/L  should  still  be  in  the  area  of  unity. 


APPENDIX  A 


Planar  Kinematic  Program  Lia  ting 

The  following  program  was  developed  to  determine  the 
equilibrium  position  of  a  four  bar  mechanism  with  a  mass 
attached  to  one  link.  The  program  calculates  the  wire 

lengths  (r2  &  r 4 )  required  to  make  the  output  link  (r3) 

level  Initially.  This  requires  the  user  to  supply  the 
lengths  of  the  input  (rl)  and  the  output  links,  an  initial 

rotation  (wint) ,  a  mlsscenter ing  of  the  optics  on  link  3  (u) 

and  the  height  of  the  mechanism  (h) ,  all  of  which  determine 
the  lengths  of  the  wires.  The  Instant  center  1/3  is  then 
positioned  directly  above  the  center  of  mass  on  the  output 
link.  A  non-linear  equation  solver  is  then  called  to  solve 
the  eight  constraint  equations  describing  the  location  ot 
the  pin  joints,  points  A  and  B.  The  equation  solver  is  a 
routine  in  the  International  Mathematics  &  Statisics  Library 
and  it  calls  an  external  subroutine  fen.  The  program  was 
written  in  Fortran  F77. 

program  maln(lnput, output, tape5»input, tape6>output) 
c 

c  Planar  model  using  Newtonian  static  equilibrium, 
c  Length  of  wires,  links  2  and  4  are  calculated 
c  based  on  an  initial  rotation  of  link  1  and  a 
c  milcentering  of  the  eg  on  link  3. 

c  Up-dated  on  Jan  2  1984 


dimension  par(l),wk(600),xo(S),xa(8) 

common/data/ rl,r2,r3,r4,r5,r6,cx,cy,dx,dy 

real  gx,gy 

external  fen, fun 

n-8 

m-  150 

pi«acos (-1.0) 
deg-180 . /pi 

Variables  describing  link  lengths  are  read 
read(5  ,*)  r3 , r 1 , di s t , ph , win t , u 
w in«wlnt/ deg 

Translate  link  1  until  the  IC  is  above  the  CG 
ax»*r3/ 2 ,+u 
ay-0 

bx-r3 / 2 .+u 
by-0 

cpx-ax*rl*(ph*cos(win)+bx*sin(wln))/(bx*(ph+rl*sin(wln))-ph*ax) 

cpy-ph 

dpx-cpx+r 1*  cos (win) 
dpy-cpy+rl*sin(wln) 
gx-0 .0 
gy-ay-dis t 

Determine  the  lengths  of  links  2,  4,  5  &  6 
r2-sqrt((ax-cpx)**2.+(ay-cpy)**2.) 
r4»sqrt((bx-dpx)**2.+(by-dpy)**2.) 
r5-sqrt((gx-ax)**2.+(gy-ay)**2.) 
r6«sqrt((gx-bx)**2.+(gy-by)**2.) 

Initial  estimate  of  independent  variables 
xo(  1 ) -ax 
xo ( 2 ) -ay 
xo  ( 3 )-bx 
xo ( 4 ) - by 
xo ( 5 )-gx 
xo ( 6 ) »gy 
xo ( 7 )=l/2 . 
xo(8)-l/2. 

Output  indicating  initial  configuration 
wrlte(6,2) 

2  format(lhl) 
write (6, 10)  rl 

10  f orraa t ( 6x , l 3hor 1 f ace  dia  -,f7.4) 
write (6, 20)  r2,r4 

20  f orma t ( 6x  ,  l 7 h leng th  of  wires  -,f7.4) 
write (6, 30)  r3 

30  f orma t ( 6x , 8hs tance  -,t7.4) 
write (6, 35)  dlst 

35  f ormat(6x , 13hplate  to  cg.-,f7.4) 
write (6, 36)  ph 

36  f orma t( 6x , 8hhe igh t  -,f7.4) 
wr 1 te ( 6 , 38 )  wint 

38  f orma t ( 6x , 2 lhinl t la  1  angular  error,  f 7 . 4) 
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MICROCOPY  RESOLUTION  TEST  CHART 
mtiowl  «um»u  of  st»no*«» 4 


wrlte(6,39)  u 

39  £ormat(6x,9hcg.  error, £7. 4) 
vrlte(6 ,40) 

40  £ormat(/// ,8x,7htheta  g  ,  6x  ,  7  A  tfte  ta  m ,  6x  ,  7  h  the  ta  2 
!  6x,7htheta  3,5x,9bd  theta  a) 

e  read  thd  number  o£  input  states  to  be 
c  solved,  1.  Read  the  input  state  tint(deg) 
read (5,5)  1 
S  £ormat( 12 ) 
do  100  j-1,1 
read(5,*)  tint 
tin- tlnt/deg 
cx-cpx+rl*cos ( tin)/ 2 . 
cy-cpy-rl*sin( tin)/2 . 
dx-dpx-rl*cos(tin)/2. 
dy-dpy+rl*sin( tin)/ 2 . 
c  Non-linear  equation  solver 

call  sscnt(fcn,lO,n,m,par,xo,£norm,wk,ler) 
c  Calculate  the  rotation  o£  link  3 
c  from  points  ay  and  by. 
ay-xo( 2) 
by-xo(4) 

tho-as in( (by-ay)/r3) 
tout-tho*deg*3600. 

Calculate  the  input  output  ratio  rat 
1  £ ( tint .ne .0 .0 )  then 
rat-tho*deg/tlnt 
else 
rat-0 
endif 

write (6, 60)  tint, tout, re t,xo(7) ,xo(8) 

>0  f orma t ( / ,  3x  ,  5  £  13 . 5 ) 

LOO  continue 

s  top 
end 


Subroutine  containing  constraint  equations 
subroutine  fen  (x,£,n,par) 

coomon/data/rl ,r2,r3,r4,rS,r6,cx,cy,dx,dy 

dimension  x(n ) , f ( n) , par ( 1 ) 

real  gx,gy 

ax*x( l) 

ay-x(2) 

bx-x (3) 

by-x(4) 

gx-x( 5) 

gy-x(6) 

£  ca«x( 7) 

£ d  b-x( 8) 


»*.  *  a  •*,  > 


•  w 


£ ( l)«(*x-bx)**2 .  +  (»y-by)**2.  -  r3**2. 

£ (2)-<ex-cx)**2.  +  («y-cy)**2.  -  r2**2. 
£<3)-<bx-dx)**2.  +  (by-dy)**2 .  -  r4**2. 

£ (4)* (*x-gx)**2 .  +  («y-gy)**2.  -  r5**2. 

£( 5)“( bx-gx)**2 .  +  (by-gy)**2.  -  r6**2. 
£(6)*(cx-ex)*£c«/r2  +  (dx-bx)*£db/r4 
£ ( 7)-(cy-ey)*£c»/r2  +  (dy-by)*f db/r4  -1. 
f (8)“( ( bx-«x)*(dy-oy)  -  (by-*y)*(dx-bx))*f db/r4 
-<gx-*x) 
re  turn 


APPENDIX  B 


Spatial  Kinematic  Program 

Tula  program  was  developed  to  perform  a  spatial 
kinematic  position  solution  for  a  five  bar  mechanism  with 
spherical  joints  and  a  mass  fixed  to  one  link.  The  input  to 
the  program  are  the  lengths  of  the  input  links  (rl,  r2 ,  r3), 
the  height  of  the  mechanism  (ph),  the  lengths  of  the  output 
links  (r4,  r5 ,  r6),  the  mlscenterlng  of  the  eg  along  the  X 
axis  (u)  and  the  vertical  distance  from  the  output  plate  to 
the  eg  (dist).  The  program  calculates  the  wire  lengths 
based  on  the  above  variables  and  an  initial  rotation  of  the 
top  plate. 

A  subroutine  (rotate)  calculates  the  final  position  of 
the  top  spherical  joints  via  a  rotation  matrix.  The  output 
rotations  are  calculated  based  on  Rodrigue's  formula  and  the 
normal  to  points  d,  e,  and  f.  As  before,  a  non-linear  equa¬ 
tion  solver  is  used  to  find  the  equilibrium  of  the  mechanism 
and  requires  subroutine  fen. 

program  main(input,output,tape5«lnput,tape6aoutput) 
real  nx  ,ny ,nz , h, bp , knx , kny , knz , kn 

dimension  xo ( 15 ) , wk( 1 104 ) , par ( 1 ) , xp ( 1 5 ) , h ( 3 ) , hp ( 3 ) 


V  v  V  V  \  V 


W-.VvWV-. 
.  *  •  -  *  *  •  »  *  ' 


* '  -* 


o  o 


i  ,ztp(3),ytp(3),xtp(3),xt(3),yt(3),zt(3) 
cooBoa/diu/rl ,  r2  ,  r3  ,  r4  ,  r5  ,  r6  ,  r7  ,  r8  ,  r9  ,  r  10  ,  r 1 1 ,  r  12 
!  ,xt,yt,*t 
«xt«ra«l  fen 


c  xd-x(l)  xe-x(4)  xf-x(7)  xcg-x(lO)  fde-x(13) 
c  yd-x(2)  ye-x(5)  yf«x(8)  ycg-x(ll)  feb-x(14) 
c  zd«x(3)  ze-x(6)  zf«x(9)  scg-x(12)  ffc-x(iS) 


n-15 

m-150 


p 1-ACOS (-1.0) 
deg-180 ./pi 

Initial  conditions,  calculate  the  geometric  center 

of  the  top  plate  and  then  points  a,  b,  and  c. 

read ( 5 , * )  rl,r2,r3 
read( 5 ,*)  r4,r5,r6 
read(5,*)  dlst,ph,u 

phll2-acos((rl**2.+r2**2.-r3**2.)/(2*rl*r2)) 
phi 13-acos((rl**2.+r 3**2. -r2**2 . ) / ( 2*rl*r3 )  ) 
phl23-pl-phll2-phil3 

phl45-acos((r4**2.+r5**2.-r6**2.)/(2*r4*r5)) 
phl46«acoa((r4**2.+r6**2.-r5**2.)/(2*r4*r6) ) 
phl56>pl-phl45-phl46 

tp>pl/2-phl23 
rp-rl/ (2*cos(tp) ) 
t b«pl/2-phl56 
r b«r4/ ( 2*cos ( tb) ) 

xtp(l)«-rp*cos(tp) 
y tp( 1)«0.0 
ztp(l)-rp*sln(tp) 
xtp(2)-0.0 
y tp( 2 )-0 .0 
*tp(2)  — rp 
xtp(3)-rp*cos( tp) 
y tp(3)-0.0 
ztp(3)-rp*sin(tp) 

Initial  estimate  for  points  d,  e,  f  and  g,  also  forces 
xo(l)"-rb*cos(tb) 
xo( 2)»ph 
xo(3)"rb*sln(tb) 
xo(4)-0 .0 
xo ( 5)«ph 
xo(6)--rb 
xo(7)»rb*cos(tb) 
xo(8)-ph 


xo(9)-rb*»in(tb) 
xo(10)«u 
xo( ll)-ph+dls t 
xo( 12 )-0 .0 
xo(13)-0. 33333 
xo(i4)-0. 33333 
xo(15)«0. 33333 

Normal  Co  lower  plate  Initially 
H  < l)-xo( 1) 
h(2)-xo(2)+l. 
h( 3)-xo(3) 

Store  for  Initial  position 
do  11,1*1,15 
11  xp(l)*xo(i) 

Initial  rotation  ot  top  plate 
read(5,*)  wxmt,wzint 

call  rotate(xtp,ytp,ztp,xt,yt,zt,wxinc,wzint) 
do  90  1-1,3 
xtp(l)-xt( 1) 
ytp(i)-yt(i) 
ztp(i)-zt(l) 

90  continue 

Calculate  wire  lengths  and  CG  constraint  links 
r7-sqrt((xod)-xtpd))**2.+(xo(2)-ytpd))**2.  +  (xo(3)-ztpd))**2.) 
r8-sqrt((xo(4)-xtp(2))**2.+(xo(5)-ytp(2))**2.+(xo(6)-ztp(2))**2.) 
r9-sqrt((xo(7)-xtp(3))**2.+(xo(8)-ytp(3))**2.+(xo(9)-ztp(3))**2.) 
rlO-sqrt((xodO)-xo(l))**2.+(xo(ll)-xo(2))**2.  +  (xod2)-xo(3)) 

!  **2 . ) 

rll«sqrt((xo(10)-xo(4))**2 ,+(xo( ll)-xo(5))**2.+(xo(12)-xo(6)) 

!  **2.) 

rl2-sqrt((xo(10)-xo(7))**2.+(xo(ll)-xo(8))**2.+(xo(12^-xo(9)) 

!  **2.) 

Input  rotation  of  top  plate,  do  (1)  times 
write (6, 170) 
write (6, 17 5) 
read(5,10)  1 
10  format(i2) 
do  350  j-1,1 
read(5,*)  omegx,omegz 

call  rotate(xtp,ytp,ztp,xt,yt,zt,omegx, omegz ) 

Call  equation  solver 

call  zscnt  ( f cn , 9 , n ,m , par , xo , t norm , wk, 1 er ) 

Rodrigue's  formula 

xha-(xo(8)-xo(2) ) *< xo ( 6 ) -xo ( 3 ) ) - ( xo ( 5 ) -xo ( 2 ) ) * (xo ( 9 ) -xo ( 3 )  ) 

yha-(xo(4) -xo( 1) )*(xo( 9 )-xo(3))-(xo(7)-xo(i))*(xo(6)-xo(3)) 

zha- ( xo ( 7 )-xo(l))*(xo(5)-xo(2))-(xo(4)-xo(l))*(xo(8)-xo(2)) 

tl-sqrt(xha**2 ,+yha**2 ,+zha**2 . ) 

hpd)-xod)+xha/tl 

hp(2)-xo(2)+yha/tl 

hp(3)-xo(3)+zha/tl 

knx- ( xo ( 2)-xp(2))*(hp(3)-h(3) )-(np(2)-h(2))*(xo(3)-xp(3)) 
kny-( hp(l)-h(l))*(xo(3)-xp(3))-(xo(l)-xp(l))*(hp(3)-h(3)) 


knz-(xo(l)-xp(l)  )*(hp(2)-h(2)  )-(hp(  l)-h(i)  )*(xo(2)-xp(2)  ) 

kn-sqr t (knx**2 .+kny**2 .+knz** 2 . ) 

ox-knx/kn 

ay-kny/kn 

az-kaz/ka 

sl-(h(l)-xp(l))*(np(l)-xo(l))+(h(2)-xp(2))*(hp(2)-xo(2))+(h(3) 

xp<3))*(hp(3)-xo(3)) 

s2-(h<  l)-xp(  1)  )*nx+(h(2)-xp(2)  ) *ny+( h ( 3 ) -xp ( 3 ) ) *nz 

s3-(h(2)-xp(2))*nz-(h(3)-xp(3))*ny 

s4-(h(3)-xp(3)  )*ax-  ( ti(  l)-xp(  1) )  *az 

s5-(h(l)-xp(l))*ny-(h(2)-xp(2))*nx 

s6-s3**2.+s4**2 .+s5**2  . 

be  ta-acos( <  sl-(a2**2  .) ) /a6 ) 

thetax-beta*nx*deg*3600 . 

the tay-be te*ny*deg*3600 . 

thetaz-beta*nz*deg*3600. 

thdm-sqrt{ thetax**2.+thetaz**2.) 

thla*sqrt( omegx**2 .+omegz**2.) 

Calculate  Input  output  ratio  rat 
if( thin. ne .0.0)  then 
rat-thdm/(thin*3600.) 
else 

rat-0.0 

Create  new  normal  if  the  tint  one  was  superseded 
h( l)-xo(l) 
h(2)-xo(2)+l. 
h ( 3 ) -xo ( 3 ) 
do  305,  1-1,13 
'05  xp(i)-xo(i) 
endif 

write (6, 174)  thdm,omegx,omegz,rat,tnorm, the tay 
'50  continue 

ou  tpu  t 

write (6, 130)  rl,r2,r3 
write (6, 160)  r4,r5,r6,ph 
write (6 ,190)  r7,r8,r9 
wrlte(6 ,135) 

write (6, 140)  (xtp( i) , 1-1 ,3 ) 
wr 1 te ( 6 , 140 )  (y  tp  (  1 )  ,  1- 1 , 3  ) 
write (6, 140)  ( z tp ( 1 ) , 1- 1 , 3 ) 
write (6, 145) 

write( 6 , 140)  ( x t( i ) , i- 1 , 3 ) 
write (6, 140)  (y t ( i ) , i- 1 , 3 ) 
wrlte(6 ,140)  { z t< i ) , i- 1 ,3 ) 
write(6 ,165) 

write (6, 140)  xo(l) ,xo(4),xo(7),xo(10) 
write(6,140)  xo(2),xo(5),xo(8),xo(ll) 
write (6, 140)  xo(3) , xo( 6 ) , xo ( 9 ) ,xo(12) 
wrlte(6,180)  wxlnt,wzlnt 

l  3  0  format! Lhl5x , 'upper  plate  dimens  ions s " , / iOx , " rl-" , 1 10 .4 , 

!  3x,',r2-%fl0.4,3x,‘’r3-',f  10.4) 


L45  forma t(/ ,  22x,'a' ,  14x,'b' ,  14x,'c') 

L40  forma t( 10x,4fL5.5) 

L  3  5  format(/,2lx,'ap',13x,'bp',13x,'cp') 

L50  format(/ ,5x, 'input  rotation  was:  (dag. )'  ,/ lOx ,2fl0 .2) 
160  format( / ,5x ,' lower  plata  dimensions lOx ,'r4«' 

I  ,£7  .3,3x,'r5-',f7 .3,3x,'r6-',f7 . 3 , 3x , ' h-' , f 7 . 3 ) 

165  format(/21x,'d',14x,'e',14x,'f',14x,'cg') 

170  forma t( I6x ,' the ta  out', 

!  ,6x  ,'omegax' ,6x , 'omegas' , 6x , 'ra t' , 6x , 'f norm' ) 

175  format(18x,'  sec.  ' , 6  x , '  sec.  ' , 6x  , '  sec.  ' ) 

174  format(/ , lOx ,6f 12 .5) 

180  format(/ ,5x, 'initial  angular  of  1 se t' , / , lOx , It  10 . 2 ) 

190  f ormat(/ ,5x,'wlre  leng ths' , / , lOx ,' r7«' , f 7 . 3 , 3x , 

!  'r8-',f7 .3,3x,'r9-', f7 .3) 
s  top 
end 

Subroutine  used  to  rotate  cop  plate 

subroutine  rotate(xl,yi,zl,xf,yf,zf, omega x , omega z ) 
dimension  xl(3),yi(3),zi(3)»xf(3),yf(3),zf(3) 
deg*180./(acos(-1.0)) 

1 f ( ( omega x+o me gaz) .It. 0.00001)  then 
do  90  1*1,3 
xf (l)-xi(i) 
yf (l)-yl(l) 
zf (l)-zl(l) 

90  continue 
else 

omegx»omegax/deg 

ooegz-omegaz/deg 

omega-sqr t( omegx**2 .+omegz**2 . ) 
ux*onegx/ omega 
uz*omegz/omega 
do  10  j-1,3 

xf(j)-xl(j)*((ux**2.)*(i-cos( omega ) )+cos ( omega ) ) 

!  +zi(j)*(ux*uz*(l-coa( omega) ) )-y 1 ( j )* s in ( omega ) *uz 
yf(j)-xl(j)*uz*sin(omega)>zi(j)*( ux*sln( omega) ) 

!  +yl(j )*cos( omega) 

zf(j)-xl(j)*ux*uz*( l- co a (omega ) )+z l( j )*( ( uz**2 . ) 

!  *(l-cos(  omega)  )-fcos(omega))+yl(j  )  *ux*s  In  (  omega) 

10  continue 
endlf 
re  turn 
end 

Subroutine  used  by  zscnt 

subroutine  fen  (x,f,n,par) 
dimension  x(n) , f (n ) , par (  l ) 

!  ,xt(3) ,yt(3) ,zt(3) 

common/deta/rl , r 2 , r3 , r4 , r 5 , r6 , r 7 , r8 , r9 , r 10  , r 1 1 , r  12 
I  »xt ,y t, zt 


f <l)-((xt<l)-x<l))*x(13)/r7)+((xt(2)-x(4))*x<14)/r8)+((xt(3) 

~x( 7 ) )*x( !5)/r9) 

£(2)«((yt(l)-x(2))*x(13)/r7)+((yt(2)-x(5))*x(14)/r8)+((yt(3) 

-x(8))*x(15)/r9)-l.O 

f (3)-((zt(l)-x(3))*x(i3)/r7)+(<zt(2)-x(6))*x(14)/r8)+<(zt<3) 
-x(9))*x(15)/r9) 

f (4)-(((x(2)-x(ll))*(zt(l)-x(3))-(x(3)-x(L2))*<yt(l)-x(2))> 
*x(13)/r7) 

+(((x(5)-x(Il))*(zt(2)-x(6))-(x(6)-x(I2))*(yt(2)-x(5))) 

*x( 14)/r8) 

+(((x(8)-x(ll))*(zt(3)-x(9))-(x(9)-x(L2))*(yt(3)-x(8))) 

*x(i5)/r9) 

£(5)"(((x(3)-x(L2))*(xt(l)-x(l))-(x(l)-x(10))*(zt(l)-x(3))) 

*x(!3)/r7) 

+(((x(6)-xU2))*(xt(2)-x(4))-(x(4)-x(iO))*(zt(2)-x(6))) 

*x( 14)/r8) 

+(((x(9)-x(12))*(xt(3)-x(7))-(x(7)-x(10))*(zt(3)-x(9))) 

*x(l5)/r9) 

f(6).(((x(L).x(10))*(yt(l)-x(2))-(x(2)-x(ll))*(xt(l)-x(l))) 

*x(13)/r7) 

+(((x(4)-x(10))*(yt(2)-*(5))-(x(5)-x(IL))»(xt(2)-x(4))) 

*x( 14 ) / r8 ) 

+  (((x(7)-x(iO))*(yt(3)-x(8))-(x(8)-x(U))*(xt(3)-x(7))> 

*x( 15 ) / r9 ) 

f(7)«(x(l)-xt(l))**2.+(x(2)-yt(l))**2.+(x(3)-zt(l))**2.-(r7**2.) 
f (8)-(x(4)-xt(2) )**2.+(x(5)-yt(2) ) ** 2 .  +  (x( 6 ) - z t( 2) ) ** 2  ( r8** 2  .  ) 

f (9)«(x(7)-xt(3))**2.+(x(8)-yt(3))**2.+(x(9)-zt(3))**2.-(r9**2.) 

f (10)-(x(4)-x(l))**2.+(x(5)-x(2))**2.+(x(6)-x(3))**2.-(r5**2.) 
f ( ll)-(x(7 )>x( l) )**2 .+(x(8)-x(2) )**2 .+(x(9)~x(3) )**2.- (r4**2.) 
f(12)-(x(7)-x(4))**2.+(x(8)-x(5))**2.+(x(9)-x(6))**2.-(r6**2.) 

£ U3)-(x(l0)-x(l))**2.+(x(ll)-x(2))**2.+(x(12)-x(3))**2.-rl0**2. 
£(14)«(x(10)-x(4) )**2.+(x( ll)-x(5) )**2.+(x(12)-x(6) )**2 ,-rll**2. 
£(15)-(x(10)-x(7))**2.+(xUl)-x(8))**2.+(x(12)-x(9))**2.-rl2**2. 
return 


appendix  c 


Planar  Stlf f ness  Program 


Tula  program  calculates  the  equilibrium  position  of  a 
four  bar  mechanism  with  flexible  members  on  links  2  &  4  and 
a  mass  attached  to  link  3.  The  user  reads  in  the  values  for 
mg,  El,  the  lengths  of  links  1,  2,  3,  and  4  (rl,  r2  ,  r3 ,  & 
r4)  as  well  as  the  input  angle  (thti  deg.)  for  each  input 
state.  As  with  the  earlier  programs  a  subroutine  is  called 
to  solve  the  9  constraint  equations  for  the  position  solu¬ 
tion.  The  program  is  written  In  Fortran  F77. 


c  Planar  analysis  with  the  wires  modeled 
c  as  small  deflected  beams.  3/16/84 
c  The  boundary  conditions  are  pinned  at  the 
c  top  and  clamped  at  the  bottom, 
c 

c  xl-tho  x4»va  x7*mb 

c  x2»ha  x5»vb  x8«da 

c  x3»hb  x6»ma  x9«db 

c 

program  main(lnput,output,tape5Blnput,tape6*output) 
dimension  x( 9 ) , par < 1 ) , wk( 1000 ) 
real  li  ,  la , 1 b  ,  lo ,mg 
integer  i,j 

common  u,d,ei,thl,tha,thb,deg,li,la,lo,lo,mg,j 
external  f cn  , cosn , sinh 
c 

pi-acos(-l.O) 

deg»180./pi 
1 1  max*1 125 
n-9 


read(5,*)  mg,ei,d,u 


« 1-4 . 94e-4 

Initial  estimate  tot  equation  solver 
x<l)-0.0 
x(2)--0.0847 
x ( 3 ) “0 .0847 
x(4)-mg/2  . 
x(5)-mg/2  . 
x(6)-0 .0 
x(7)-0.0 
x (8 )-0 . 0 
x(9)-0.0 

User  Inputs,  do  (1)  times 
read(5,2)  1 
format( 12 ) 
do  40  1-1,1 

read(5,*)  li  ,  la , lb , lo , tntl 
thi- th  tl/deg 

Calculate  initial  angles 
1 f < 1 . ne . 1 )  goto  5 
xa-(  lo- 11 ) / 2 . 
tha-asln(xa/la) 
ttib-asln(xa/lb) 

5  continue 

Call  equation  solver 

call  zscnt( f cn , 13 ,n , 1 traax , par , x , f norm ,wk , ler) 
Calculate  output 
da«x( 8) 
d  b-x( 9) 
pba-tha-thi 
pbb-thb+thi 
tot»x( 1 ) *deg 

Calculate  Input  output  ratio 
rat-0.0 

If  (  tlitl .  le  .0 .0)  goto  11 
rat-tot/  thtl 
11  continue 

xb-ll*coa(thl)+lb*sin(phD)-dD*cos(pnb) 
yb-li*sln(thl)-lo*cos(phb)-db*sin(pbo) 
xa--la*sln(pba)-da*cos(pba) 
ya--la*cos(pha )+da*sln(pha ) 

Output  critical  variables 
vr 1 te ( 6 , 20 )  th tl , to t* 3600 . , ra t 
vr lte ( 6 , 20 )  li , lo , la , lb, el 
10  format( 10x,4el5 .4) 
vr 1 te ( 6 ,30 )  (x(j),j-l,7) 
vrl te( 6 ,30 )  (x(J),j-8,9) 

write <6, 30)  xa,ya,xb,yb,pba*deg,pbb*deg 
30  f orma t( 5x , 7f 1 2 . 5 ) 

20  f orma t( 8x ,7 f 10 .4 ) 

J-o 

40  continue 


Subroutine  used  by  mat 

subroutine  fen  (x,£,n,par) 
dimension  x(a),f{n)lpar(l) 

common  u,d,el,thi, the ,tnb,deg,li,la,lb,io,mg,j 
real  11 , La , Lb , lo ,ma ,mb ,mg , ha , ftb 
external  cosh,sinh 
integer  l,J 

tho»x( l) 
ha«x ( 2 ) 
hb*x(3) 
va«x(4) 
vb-x( 5) 
ma-x( 6 ) 
m  b*x( 7 ) 
da«x( 8) 
db*x(9) 

pha» tha- thl 
phb-chb+thi 
bl*aqrt( va/ el) 
b2«*aqr  t (  vb/  e  1 ) 
cblla«coab(ol*la) 
abLLa»alnh(bl*ia) 
cb2Lo>coah(o2*lb) 
sb2Lo«alnh(b2*lb) 

£(l)»ha*cos(pha)+va*siu(pha)+hb*cos(phb)-vb*8ln(phb) 
f  (2)»v«*coa(pha)-h»*sln(p6a)+ti6*»ta(phb)+vb*co8(pnb) 
£ ( 3  )*-ma-mo+vb* lo*cos(phb-t ho )+bb*lo*sln(phb-tho) 
-mg*((Lo-u)*cos(tho)/2 .+d*a ln( tho) ) 

f(4)»ha*abLLa/ ( bl*cbl la )+ma 
f (5)»hb*ab2lb/( b2*cb2lb)+ob 

f(6)»-ha*ablla/(va*bl*cblla)+( tho+ha/ va)*la-da 
£(7)»-no*sb21b/(vb*b2*cb21b)  +  ( tho+h  b/vb)*lb-db 

xa*-La*sin(pha)-da*cos(pha) 

ya*-la*cos(pha)+da*ain(pha) 

xb«li*cos(thl)+lo*slu(phb)-db*cos(phD) 

yb*ll*sln(tftl)-lb*cos(pftb)-db*sia(phb) 

£(8)“(xb-xa)**2.+(yb-ya)**2.-lo**2. 

f(9)*(yb-ya)/(xo-xa)-tan(tho) 

return 

end 


Evaluate  the  slab  and  cost)  in  function] 


function  cosb(a) 


APPENDIX  D 


Project  Staff  in  1983-1984 


Faculty 


M.  M.  Barash.  Bans  burg  Professor  of 

Manufacturing  ana  Professor  of  Industrial 

Engineering .  Principal  Investigator  & 

Project  Director 


C.  R.  Liu.  Professor  of  Industrial 

Engineering . 

K.  S.  Fu.  Goss  Distinguished  Professor 
of  Engineering  (Elec.Eng.) . 

J.  Modrey.  Professor  of  Mechanical 
Engineering . 

A.  L.  Sweet.  Professor  of  Industrial 
Engineering . 


Principal  Investigator 
Faculty  Associate 
Co— Principal  Investigator 
Co-Principal  Investigator 


W.  Stevenson,  Professor  of  Mechanical 
Engineering . . 


Faculty  Associate 


J.  J.  Talavage.  Professor  of  Industrial 

Engineering .  Faculty  Associate 

R.  Hannam*.  Visiting  Associate  Professor 

of  Industrial  Engineering .  Faculty  Associate 

W.  Johnson"*,  Visiting  Professor  of  Industrial 

Engineering .  Faculty  Associate 


A.  Shumsherrudint,  Visiting  Associate  Professor 

of  Industrial  Engineering .  Faculty  Associate 


(  "University  of  Manchester  Institute  of  Science  and  Technology,  England 
""University  of  Cambridge,  England 
tCranfield  Institute  of  Technology,  England) 


Graduate  Resear  oh  Assistants 

P.  Chen 
Y.C.  Chou 
P.  Ferreira 

R.  Khanna 

S. K.  Lee 
G.R.  Liang 
Y.T.  Lin 
J.  Lopez 
D.  Noller 


Y.S.  Ouyang 
U.  Roy 
S.  Shodhan 

R.  Srinivasan 

S.  Venkatramen 
R.  Venugopal 
M.C.  Wu 

J.  York 


